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I INTRODUCTION 


THE APL ECONOMETRIC PLANNING LANGUAGE is an interactive 
planning language based upon APL. It provides features for 
dealing with economic variables (primarily time-series), such 
as data analysis and transformation, tabular and graphic display, 
parameter estimation, model solution and file handling. 


The APL ECONOMETRIC PLANNING LANGUAGE provides the prac- 
ticing economist, or the teacher, with easy-to-use tools for 
interactive model building and solving. Care has been taken 
to use much of the power inherent in APL, such as its rich set 
of primitive functions and its compact data handling and logic 
capabilities. Operators on time series have been designed so 
as to be compatible with the underlying language. 


As a consequence, the APL ECONOMETRIC PLANNING LANGUAGE 
constitutes a truly interpretive and interactive system, with 
all attendant advantages. The objects with which it deals are 
APL functions and variables. Therefore, the user can easily 
add his own functions to tailor the language to his specific 
requirements. 


For heavy usage, data file routines have been added, based 
upon the processor TSIO of APLSV or AP110 of APL/CMS. The user 
can set up and utilize his own data base by means of simple store, 
retrieve and update functions. Those users who already have data 
bases, should find little difficulty in adding their own interface 
functions via TSIO or AP110. 


The system was designed and developed at the IBM Philadelphia 
Scientific Center (see refs. 10-12), with consultation from IBM's 
Business Planning Department, Stuttgart, and Economic Research De- 
partment, Armonk. It is installed at the IBM Economic Research De- 
partment, Armonk, and the APL Design Group, Palo Alto. 


The user of the APL ECONOMETRIC PLANNING LANGUAGE may wish 
to become somewhat familiar with APL itself. Relevant documenta- 
tion on APL is available in the APL/360 Primer (GH20-0689), the 
APL/360 User's Manual (GH20-0906), and the APL Language Manual 
(GC26-3847). 


II SYSTEM OVERVIEW ወ 


The system comes in the form of two APL workspaces on 
one tape: ΕΡΤ, EPL2. It was developed and can be used in 
one 80K byte workspace, but has been regrouped into two work- 
spaces of size less than 65K. 


The functions of the system are partitioned into seven 
groups. 


EPL1 contains: 


BASEGRP 
DISPLAYGRP 
PLOTGRP 


EPL2 contains: 


FILEGRP - for use with APLSV 
FILEGRPC - for use with APL/CMS 
REGRESSGRP 

SOLVEGRP 

TRANSLATEGRP 

ADDENDUM 


The user will work with an active workspace which will 
contain all those groups which he needs at a particular time. 
The active workspace must always contain BASEGRP. The other 
groups are "copied" (see III OPERATIONS) as needed. They con- 
tain those functions which are available to the user: 





DISPLAYGRP: DISPLAY TABULATE 
PLOTGRP: PLOT PLOTF VS 
FILEGRP: CDATA DDATA EDATA RDATA RETRIEVE STORE 
FILEGRPC: CDATAC DDATAC EDATAC RDATAC RETRIEVEC STOREC 
REGRESSGRP: 44 COR DIST GLS INST PCINST REGRESS 
SOLVEGRP: ORDER SOLVE SQLVECON SOLVEIT SOLVEREL 

SSOLVE SSOLVEF STRUCTURE 
TRANSLATEGRP: DO TRANSLATE 
ADDENDUM: RETRANS VOCABULARY UPDATE PREP TABULTEXT 


In addition, these groups may contain auxiliary functions 
which can not be called directly by the user and are not described 


in this report. o 























The user will normally work on a 2741 terminal with 
APL type ball 987, and ህ[11 have no difficulty in composing 
function names with understruck capital letters. All 
functions listed above are described in this manual to an 
extent sufficient for their use. If the user at some point 
wants to make changes and additions to the system, he may 
wish to consult the SYSTEMS GUIDE. 


In the systems guide all functions are listed , most 
of them alphabetically. The alphabetic order is not observed 
for the functions described in the ADDENDUM of the PDOM (the 
functions in the group ADDENDUM of EPL2). They were added as 
an afterthought and are listed together at the end of the 
systems guide. 


In a sense this is not undesirable. The addition of 
these functions provides the user with an example of the 
extenstbtltty of the system. 
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3.1 


PROGRAM INSTALLATION 


For use under APLSV, the APL Econometric Planning Language 
workspaces are loaded into the APLSV public library using 
the system-supplied utility. Under VM/370 CMS, the pro- 
gram's APLSV workspaces can be converted to APL/CMS work- 
spaces and loaded into the VM/370 system, using the system- 
supplied utility. For use under VS APL, the APL/CMS work- 
Spaces are converted to VSAPL using the system-supplied 
utility. 


The file-handling group FILEGRP will operate under APLSV; 
the file-handling group FILEGRPC will operate under CMS, in 
either the APL/CMS version or the VS APL version running 
under VM/370 CMS, but not in VS APL under a VSPC environ- 
ment. Except for file-handling, the program will work in 
VSAPL under a VSPC environment. 

















LLL OPERATION 


The system has been designed such that very little 


knowledge of APL ts required for its operation. The user who 
has no famillarlity with APL can start a terminal session in 
the followIng fashion: | 


1. 
2. 


Load the workspace EPL1 by the command: )LOAD EPL1 
Change the (active) workspace name to, say, TEST: 
)WSID TEST. 


Work In the active workspace TEST by following the 
examples In thls manual, If desirable, save the active 
workspace by the command: )SAVE. | 


"copy" other functions from EPL2 as required, ege: 


)COPY EPL2 REGRESSGRP to obtain all functions of the 
group REGRESSGRP, or | 


)COPY EPL2 INST to obtain only the function INST of the 
REGRESSGRP, 


Groups or functlons (or variables) no longer needed can 
be removed from the active workspace by commands such as: 
)ERASE DISPLAYGRP or )ERASE TABULATE. 





IV SYSTEM DESCRIPTION | @ 


This manual describes the APL ECONOMETRIC PLANNING 
LANGUAGE. It is assumed that the reader is familiar with the 
principles of econometrics, but references to the relevant 
literature will be given in the text. The APL ECONOMETRIC 
PLANNING LANGUAGE concentrates on parameter estimation for 
multivariate models and on numerical solution of these 
models. According to the usual steps in practical 
econometric work the report is divided into the chapters 


TIME SERIES: DEFINITION AND DISPLAY 

MODEL LANGUAGE 

REGRESSION 

MODEL SOLUTION 

DATA FILE ROUTINES 

To demonstrate the capabilities of the APL ECONOMETRIC 

PLANNING LANGUAGE an example is repeatedly used. Context and 
data are taken from BASS [2], but the formulation of the 
model is changed. The model describes the interrelation between 
Sales and advertising expenditures in the US Cigarette industry. o 
SAMPLE MODEL 


a(1,2)51(t-1) + 


(1) 51(1) = a(1,1)W1(t)/W2(t) 4 
a(1,3)DI(t) + u(1,t) 

(2) S2(t) = a(2,1)W1(t)/W2(t) + a(2,2)S2(t-1) + 
a(2,3)DI(t) + u(2,t) ` 

(3) Wl(t) = a(3,0) + a(3,1)S1(t)/S2(t) 4 
a(3,2)S1(t-1)/S2(t-1) + u(3,t) 

(4) W2(t) = a(4,0) + a(4,1)S1(t)/S2(t) + 
a(4,2)S1(t-1)/S2(t-1) * u(4,t) 


sales per capita of filter [nonfilter] 
cigarettes 

years 1953 to 1965 

advertising expenditures per capita for 
filter Cnonfilter] cigarettes ($) 


S1(t), rS2(t)j 


W1(t), r[W2(t) ገ 





DI(t) 
u(n, t) 


disposable personal income per capita ($) 
disturbances. 


W1Ct), W2(t), DICt) are deflated ኩሃ appropriate 
price indices. The population comprises all persons over 20 


years of age. 


TABLE 1: Data for the sample model (BASS). 


Year 51 $2 W1 W2 DI 

1953 250.228 3196.318 «055 .520 2609.336 
1954 398,658 2856.406 «126 .425 2622.769 
1955 690.081 2668.824 .236 «271 2785.544 
1956 952,525 2453.522 .367 «260 2895.810 
1957 1939. 111 2178.211 .453 .279 2914.137 
1958 1414.718 2008.075 .451 .289 2904.290 
1959 1526.195 2006.827 .495 ,287 3018.978 
1960 1571.303 1991.590 .462 .292 3055.062 
1961 1613.578 1971.696 .455 .286 3115.875 
1962 1657.4837 1900.684 .434 .283 3221.365 
1963 1731.530 1808.007 «518 .285 3293.679 
1964 1672.207 1590.523 ,506 .283 3472.402 
1965 1702.119 1633.202 .486 «217 3658.895 








The technical part of econometric modeling consists of 
1. the statistical estimation of the unknown values for the 
parameters a(1,1)...a(4,2) and 2. the solution of the model 
for the endogenous variables S1(t), S2(t), W1(t) and W2(t), 
when values for the parameters and the exogenous variables 
DI(t), S1(t-1) and S2(t-1) are provided. A solution can be 
computed either for historical values of the exogenous 
variables (model validation) or for estimated future values 
(forecasting). | 


The system provides limited information techniques for 
the estimation of models that are linear in the parameters, 
i.e. ordinary regression, generalized regression, 
instrumental variables substitution with or without 
principal components transformation, and polynomial 
distributed lags. Models may be nonlinear in their 
variables, as is the sample model. Little emphasis is put on 
the analysis of the internal structure of a single time 
series. The system is modular so that either the user's own 
APL functions or public library functions for time series 
analysis can be added easily. | 





1. TIME SERIES: DEFINITION AND DISPLAY 


ED CD Oe 0 eee 


In a general sense, data in econometric modeling may 
be time series, cross-sectTonal observations, equations or 
complete models. In this chapter we will consider the 
definition of time series (and cross-sectional series) and 
various display routines. From an APL potnt of view time 
series will be represented as numeric arrays, while 
equations and models will have the form of character arrays; 
see chapter 2. Optional names can be assigned to each data 


object. These names should not contain understruck 
characters, 


1.1. DEFINITION OF A TIME SERIES 


In discrete econometric modeling equidistant time 
series are used almost exclusively. Therefore, it is 
sufficient to concentrate on this type of series. A time 
series with an optional NAME is defined by the dyadic 
function DE 


NAME*HDR DE VAL 


The right argument VAL contains the values of the 
different periods, ordered with increasing time. VAL ts 
usually a vector, denoting a series with exactly one value 
for each period. But VAL can also be a higher rank input, 
denoting a series with multiple classification. The last 
coordinate must always relate to time . As an example 
consider a compound series SALES, where the first coordinate 
relates to different products, the second to geographically 
distinguished markets and the third to different sales 
periods. 


The time characteristic of a series is described by 
the left argument HDR, a three-element vector, subsequently 
called header. HDR[1] is the periodicity of the series and 
is defined as the number of periods corresponding to one 
full year. The APL ECONOMETRIC PLANNING LANGUAGE recognizes 
the following valid periodicities: 


1 = yearly data 

2 = semi-yearly data 
4 = quarterly data 
12 = monthly data 
52 = weekly data, 











The other two elements of the left argument HDR 
determine the origin of the series, which is the first 
period of the VAL array. HDR[2] is the year and HDR[3] the 
period within the year. The following are valid examples of 
headers: 


4 1970 2 quarterly series starting 2. Q. 1970 
12 1970 5 monthly series starting May 1970 
1 1971 1 yearly series yearly starting 1971 


An error message DF SPEC ERR is returned if the header 
is not correctly specified, 


The explicit result WAME is a global APL-variable and 
can be manipulated outside the APL ECONOMETRIC PLANNING 
LANGUAGE by any APL-primitive or function, NAME is basically 
the catenation of the header and the value array. For 
internal system conventence the header is decoded into a 
scalar. If needed, the user can encode it again by applying 
the function EC? see example below. 


The function 
LIST ASSIGN X 


assigns names to the rovs of a rank two series X (matrix), 
so that an aggregate series can be decomposed into several 
rank one series (vectors). LIST contains those names. An 
ASSIGN SPEC ERR message is returned if the number of rows in 
X does not correspond with the number of names in LIST. 
After the execution of ASSIGN, the various rows can be 
called and manipulated individually by their new names. 


Example 1 
1.2. CROSS-SECTIONAL SERIES 


For dynamic econometric models, repetitivity of 
economic phenomena is associated with the time component. 
There are, however, also time invariant models with 
observation data sampled in a cross-sectional manner. This 
data type can also be defined by means of the function DF 
with a scalar zero header H#DR<+0. In the following we will 
omit the treatment of cross-sectional series, since the main 
focus of the APL ECONOMETRIC PLANNING LANGUAGE is on dynamic 
modeling. But most of the functions work also for 
cross-sectional data (where this is meaningful). 





 COST*«12 1973 5 DF 4.1 5.2 5.9 6.3 7.1 6.8 7 5.8 5.1 


COST | 
2368112 4.1 5.2 5.9 6.3 7.1 6.8 7 5.8 5.1 


ጀሮ COST[1] 
12 1973 5 


SALHS«4 1972 3 DF 2 2 362 4 3 7 5 6 3 2 11 9 13 8 





SALES 
789104 2 4 3 
789104 7 S 6 
789104 3 2 11 
789104 9 13 8 


TEST+1 1965 1 DF 3 497 3 u 9 11 9 5 7 3 17 15 13 


TEST 
196601 7 3 4 9 
196601 11 9 5 7 
196601 3 17 15 13 


'TEST1,TEST2,TEST3" ASSIGN TEST 


TEST 1 
196601 7 3 4 9 


TEST 2 
196601 1195 7 
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1.3. TABULAR DISPLAY OF TIME SERIES 


The monadic function 


DISPLAY LIST 


displays all time series contained tn the right argument 
LIST, As with all defined. which 
accept a LIST of series as. arguments, these series must be 
of rank one (vectors), The only exceptions are the functions 
in chapter 5, DATA FILE ROUTINES. 


Example. 2 


A simultaneous tabular display of several time series 


» 1s provided by the function 


CTR TABULATE LIST 


The right argument LIST is a list of rank one series, 
which must.all have the same periodicity. The tabulatton has 
as domaln those periods common to all series (intersection 
of the series; for explicit contro] of the intersection 
process by the user; see chapter 1.5.). Blank lines between 
successive series are inserted by additional commas in LIST; 
see example. 


An INTERSECT PER ERR message Is given [የ 
periodicities do not match. No output is returned if the 
intersection of the series is empty. 


The left argument CTR ts a h- -component vector of 
parameters: CTR[1] ts an indent in front of the first listed 
period. It should be at least as large as the longest name 
in LIST. CTR[2] defines the number of periods within a table 
page. Set CTR[2]<0 if the number of periods per page should 
be calculated automatically. ZABULATE then assigns as many 
periods to one page as will fit according to the maximally 
permitted printing width (OPW), CTRL3] and CTRL^] control 
the format of the tabulated values. They should be chosen so 
that all significant positions will be printed. CTR[3] is 
the total number of positions, CTR[4] the number of 
positions after the decimal point. If the absolute value of 
a table entry is too large, a blank is printed instead. 


A TABULATE SPEC ERR message is returned if CTR is not 


correctly specified. 
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DISPLAY 'COST,S1' 


VARIABLE COST 


PERIODICITY. = 12 
ORIGIN = 1973 5 
NO OF ENTRIES = 9 





TIME VALUE TIME VALUE 
1973 5 4.10000 1973 10 6.80000 
1973 6 5.20000 1973 11 7.00000 
1973 7 5.90000 1973 12 5.80000 
1973 8 6.30000 1974 1 5.10000 
1973 9 7.10000 
VARIABLE S1 
PERIODICITY = 1 
ORIGIN - 1953 1 
NO OF ENTRIES = 13 
TIME VALUE TIME VALUE 
1953 1 250.32800 1960 1 1571.30300 
1954 1 398.65800 1961. 1 1613.57800 
1955 1 690.08100 1962 1 1657.48700 
1956 1 952.42300 1963 1 1731.53000 
1957 1 1232.11100 1965 1 1672.20700 
1958 1 1414.71800 1965 1 1702.11900 
1959 1 1526.19500 
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EXAMPLE 3 


10 5 9 3 TABULATE 'S1,S2,W1,W2,,DI' 











1953 1951 1955 1956 1957 
Si 250.328 398.658 690.081 952.423 1232.111 
22 3196.318 2856.406 2668.824 2453.522 2178.211 
W1 .055 .126 .236 .367 453 
W2 .520 2425 .371 .360 .279 
DI 2609.336 2622.769 2785.544 2893.810 2914.137 
1958 1959 1960 1961 1962 
S1 1414.718 1526.195 1571.303 1613.578 1657.487 
S2 2008.075 2006.827 1991.590 1971.696 1900.684 
W1 451 495 462 455 434 
W2 .289 .287 . 292 .286 .283 
DI 2904.290 3018.978 3055.062 3115.875 3221.365 
1963 1964 1965 
S1 1731.530 1672.207 1702.119 
52 1808.007 1590.523 1633.202 
W1 .518 .506 .186 
W2 .285 .283 . 317 
DI 3298.679 3472.402 3658.895 
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1,4, PLOTS OF TIME SERIES 


The plot functions incorporate several functions from 
Graphs and Histograms In APL [5]. | 


The function 


CTR PLOT LIST 


plots all series contained in LIST within the same graph. 
The abscissa is time. If instead of time, a time series x [5 
to be used as abscissa, it must be indicated by the use of 


፻2፡ 
CTR PLOT LIST ፻፪ ιχὶ 


An INTERSECT PER ERR message Ís produced [የ the | 
periodicities of the series tn ¿IST do not match. Also, χ 
must match if vs is applied. Again as in the previous 
chapter only the common Intersection of the series is 
plotted. | | 


The three element vector CTR controls the size and 
location of the graph. CTRL1] is the length of the ordinate 
(In printed vertical lines), CTRL21] the length of the 
abscissa (in printed horizontal positions) and CTRL3] an 
‘indent to the ordinate. If the user wants to produce a finer 
plot with typing element #1167114, the function PLOTF should 
be applied instead of PLOT. However, in this case LIST can 
contain a maximum of only two series. A PLOTF SPEC ERR 
message Is returned if LIST refers to more than two series. 
Note also that the fine plet can be used only on a 2741 with 
correspondence keyboard. u 


Example 4 


1.5. EXPLICIT INTERSECTION CONTROL 


In the above examples of TABULATE and PLOTF the result 
extended over all common peridos of the specified list of 
time series (largest common intersection). The APL 
ECONOMETRIC PLANNING LANGUAGE provides the user with an 
explicit control capability for the intersection process. He 
can specify the time frame to which the intersection process 
1s to be restricted. However, if the specified time frame 
includes the largest common 
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@ EXAMPLE 4 


50 40 10 PLOTF 'S1,52' 




















ABSCISSA = TIME STARTING FROM 1953 
. = S1 
+ = 52 
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intersection, the latter ts taken automatically. The time | @ 
frame capability allows the user to limit his analysis to a 

selected subset of the avallable data. The advantage of such 

a feature wlll become obvious In the later sections 

REGRESSION and MODEL SOLUTION. 


A time frame is set with the dyadic function ST 
HDR1 ST HDR2 | 


HDR1 and HDR2 are the first and last periods (inclusively) 
of the frame. The arguments are three- element headers, as 
described in chapter 1.1. A specification HDR1+*0 sets no 
explicit lower bound, but imposes that this bound be defined 
by the largest common intersectiom. Similarly, HDR2-0 sets 
no explicit upper bound. The headers can be specified in any 
periodicity independent of the periodicity of the series 
under consideration. Internally equivalent bounds are 
computed for all possible periodicities., | 


Once set, a time frame remains active for all 
subsequent operations. It can be erased by 


ο ST 0 


From then on the largest common intersectlon Is taken until] @ 
a new frame is specified. A DF SPEC ERR message occurs [የ 
one of the headers HDR1 or HDR2 ts misspelled. 


The current frame setting can be queried by execution 
of the niladic function 


FRAME 


It returns the lower and upper bounds for each perlodicity; 
or a zero if no frame has been set. 


The Intersection of two time series X and Y ts agaln a 


(possibly shorter) time series. It can be found with the 
function 


(NAME- ) X IS Y 
X and Y must have the same rank and dimension with the 


exception of the last coordinate (time). The rank of the 
result is Increased by one. 





The tntersection of a LIST of rank one series ts 
obtained by 
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( NAME+ ) INTERSECT LIST 


Both IS and INTERSECT return an empty vector if an 
error [5 detected or the Intersection fs empty. In the 
latter case IS causes the message IS EMPTY ERR. IS delivers 
an IS RANK ERR tf X and Y do not have the same rank. An 
INTERSECT SPEC ERR for INTERSECT Indicates that LIST 
contains no non-empty series. For both IS and INTERSECT an 
INTERSECT PER ERR ts produced if the periodicities of the 
serles do not match. 
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EXAMPLE 9 


1 1958 1 


IN 
|-3 


1 1961 1 


10 0 9 3 TABULATE 'S1,DI' 


1958 1959 1960 1961 
23 1414.718 1526.195 1571.303 1613.578 
DI 2904.290 3018.978 3055.062 3115.875 


O ST 12 1956 4 


10 0 9 3 TABULATE 'S1,DI' 


1953 1954 1955 


ጋ 250.328 398,658 690, 083 
DI 2609.336 2622.769 2785.544 


12 1956 4 
52 1956 17 


O O O O O 
O O O O O 
O O O O O 

J- 

= 

ርጋ 

on 

Ο 

= 


FRAME 
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1 1963 1 ST O 


S+S1 IS S2 

S 
196501 1731.53 1672.207 1702.119 
196401 1808.007 1590.523 1633.202 


DATA+INTERSECT 'S1,52,W1,W2,DI' 





DATA 
1.96450100E5 1.7315300E3 1.6722070E3 1.7021190£23 
1.9640100F5 1.8080070E3 1.5905230£3 1.633202053 
1.9640100£5 5.1800000F 1 5.0600000F 1 4u.,8600000F 1 
1.96401002F5 2.8500000F 1 2.8300000F 1 3.170000072 1 


1.9640100£5 3,2986790ኛ3 3.4724020E3 3.658895073 
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2. MODEL LANGUAGE 


The APL ECONOMETRIC PLANNING LANGUAGE has several 
arithmetic and logical operators which In effect make most 
of the APL primitive operators available for the 
manipulation of time series, i.e. arrays associated with a 
header. It also has some special operators frequently needed 
in econometric modeling. This permits the user to formulate 
compound expressions and models. Due to the special 
structure of the operands, APL operator symbols cannot be 
used directly but must be replaced by mnemonic functions 
such as P for + or M for -, etc. However, for the user who 
prefers symbols, a translation routine is provided to 
generate the mnemontc function names before further 
execution; see 2.5. 


Two situations can be distinguished, Efther an 
operation on one or several time series is carried out 
immediately in order to define a new series (e.g., taking 
the logarithm of a raw series and subsequently using only 
this transformed series), or a model is to be formulated for 
repeated execution. In both cases the same operators are 
used, but in the latter case the operations have to be 
stored, In the following paragraphs 2.1.,2.2. and 2.5. we 
shall discuss the definition and the immediate execution of 
operations, in 2.4. we shall treat the formulation of 
models. The execution of models will be partly covered in 
2.5. (namely for a specific single equation case), but more 
extensively in chapter h. 


2.1. BASIC OPERATORS OH TIME SERIES 


A dyadic operator with a symbolic name F 
( ΝΑΗΕ- ) X F Y 
returns the following result: 


a) X and Y are both time series. The operation F will be 
performed for the last coordínate (time) of the 
intersection x ZS Y (see 1.5.). A possible setting of a 
time frame is observed. An optional NAME can be assigned 
to the result to define a new series, or the result can 
be immediately used as Input to a subsequent operation; 
see examples below. Perlodicity and rank of X and Y must 
be equal, also the dimension vector except for the last 
coordinate, Otherwise intersection error messages will be 
printed and an empty vector result returned; see 1.5. 
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b) Either X or Y is a time series and the other argument a 
scalar. The operation takes place between all values of 
the series and the scalar. A time frame setting is not 
observed, | | 

ο) Both X and Y are scalars. The operation is performed 
between the two scalars in the usual way. This 
"degenerate case" might be replaced by the direct use of 
the APL primitive symbol. 


A monadic operator with the symbolic name H 
( NAME- ) E X 


accepts either a time series or a scalar argument X. In the 
first case the operation [5 executed for all values of the 
series. A time frame setting is not observed. To a scalar X 
the operator is applied in the usual way and again might be 
replaced by the corresponding APL primitive symbol. 


Several operations can be combined within one 
expression. The evaluation is done strictly from right to 
left, if not differently speclfled by the use of 
parentheses. 


Table 2 contains a summary of all operators provided 
by the APL ECONOMETRIC PLANNING LANGUAGE, 


Example 6 
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EXAMPLE 6 
0 ST 1 1956 1 


Si P S2 
195401 3446.646 3255.064 3358.905 3405.945 
WL-LOG Wi D #2 
WL - - 
195401 ~2.2464956 “1, 2158073 0.45237026 0.019257817 
LL-0 GE LOG W1 D W2 ο 


LL 
195401 1 1 1 O 


ይይሩ0 GE LOG #3 D #2 


LL 
195401 1110000000000 
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TABLE 2: Summary of operators defined in the APL ECONOMETRIC 
PLANNING LANGUAGE 





System 
Operation Function 
Name 
Diadic + P 
S M 
x I 
፥ D 
* PW 
< L 
< LE 
S E 
2 GE 
d G 
z NE 
A AND 
v OR 
@ Monadic log LQG 
DI exp EXP 
sin SIN 
cos COS 
tan TAN 
abs value ABS 
cumul value CUM 
I 5t NOT 
Spec. Oper. shifts LAG 
x first diff. DEL 
first quot. RTQ 
rect.distr. BECT 
norm.distr. NORM 


change period. CHANGE 





23 


2.2. SPECIAL OPERATORS 


There are some special operators which are frequently 
used in the formulatton of econometric models. 


Shift operator: Shifts in time are obtained by the 
application of the ZAG functton 


( NAME- ) K LAG X 


X ts the sertes to be shifted. A positive integer X produces 
a lag of K periods, a negative integer X a lead of X 
periods. This monadic operator does not observe a time frame 
setting. 


Elrst Differences: The function 
( NAME< ) DEL X 
returns the first difference of X which is defined as 
X M 1 LAG X 


Since this is essentially a dyadic operator (M), DEL does 
observe a time frame setting. 


First Quotients : The function 
( WAME+ ) RTQ X 
returns a result equivalent to 


X D 1 LAG X 


A time frame setting is observed, 


Rectangular Random Numbers: The function 
( NAME+ ) RECT K 


evaluates K random numbers, untformly distributed within the 
unit interval (0,1). If a time series with random number 
values is to be defined, an approprtate header must be 
catenated (for DF see 1.1.) : 


NAME*HDR DF RECT K 


Normal Random Numbers: K random numbers with standard normal 
distribution N(0,1) are obtained by 
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( WAME+ ) NORM K 


Changing the Periodicity of a Lime Series: If it becomes 
necessary to generate, e.g., quarterly series from monthly 
serles or vice versa, or tn general to change the 
perlodicity of a series, this can be done by means of the 
CHANGE function: 


( NAME- ) CTR CHANGE X 


CTR ås a two-element control vector. CTRL1] denotes the 
perlodicity of the new series (NAME), If tts pertodicity is 
smaller than that of the old series X, several periods of X 
have to be compressed into one new period. By CTRL[2], one of 
the following three compression techniques can be selected: 


CTRL2] *1 average value 
2 cumulative value 
3 1851 period value. 


| If the new periodicity Is larger, i.e. the series 
should be expanded, a choice between ten different 
Interpolation techniques can be made: 
CTRL2]- 1 stepwise interpolation 
2 linear interpolation 
3 interpol. by compound growth rate 
4 equal parts 
S sum of digits decreasing 
6 sum of digits increasing 
7 zeros 
8 cumulative equal parts 
9 cum. sum of digits decr. 
10 cum. sum of digits incr. 

Techniques 1, 2 and 3 interpolate between two 
successive values of the old series; techniques H, 5 and 6 
spread the values over the inserted periods; technique 7 
sets zeros Into the inserted periods; techniques 8, 9 and 10 
spread according to h, 5 and 6, and cumulate. 


| El ther the old (compression) or the new (expansion) 
periodicity must be an integral multiple of the other 
periodicity . Otherwise a CHANGE PER ERR message is 
delivered. CHANGE SPEC ERR indicates a misspelling of 
CTRUL2]. A CHANGE INT ERR occurs if interpolation technique 2 
or 3 is attempted for a single-period time series X. 


Example 7 
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EXAMPLE 7 


0 ST 1 1956 1 


ES1-(227.2 T W1 D W2) P (.55 


IN 
pa 
ኩ 
I> 
Ic 


ES1 
195501 519.7708 698.05333 958.41953 


RANDOM-u 1970 3 DF NORM 7 


RANDOM _ _ | 
788304 0,5037507 0.42456014  2,.0418944% 0.30000154% “0.560151 99 
0.92051797 0.33819378 


Γκι 1970 3 DF 15 20 30 50 100 70 40 10 


T11+1 1 CHANGE T 
T12-1 2 CHANGE T 
T13-1 3 CHANGE T 
T21-12 1 CHANGE T 
T22-12 2 CHANGE T 
T23-12 3 CHANGE T 
T2u-12 4 CHANGE T 
T25+12 5 CHANGE T 
T26+12 6 CHANGE T 
T27-12 7 CHANGE T 
T28<12 8 CHANGE T 
729-12 9 CHANGE T 
7210-12 10 CHANGE T 
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10 0 8 2 TABULATE 'T11,712,713' 


1970 1971 1972 
T11 17.50 62.50 25.00 
T12 35.00 250.00 90.00 
T13 20.00 70.00 10.00 


0 ST 1 1970 1 


LIST+'T21,722,723, ,724,725,726, ,727, ,728,729,7210! 











T21 
T22 
129 


T24 
125 
126 


T27 
T28 
T29 
T210 


1970 7 


15.00 
11.67 
12.13 


9.00 
7.50 
2.50 


00 
5.00 


7.50 
2.50 


10 6 8 2 TABULATE LIST 


1970 8 


15.00 
13.33 
13.49 


5.00 
5.00 
5.00 


. 00 
10.00 


12.50 
7.50 


1970 9 


15.00 
15.00 
15.00 


5.00 
2.50 
7.50 


15.00 
15.00 


15.00 
15.00 


20.00 
16.67 
16.51 


6.67 
10.00 
3.33 


. 00 
6.67 


10.00 
3.33 


20.00 
18.33 
18.17 


6.67 
6.67 
6.67 


«00 
13.33 


16.67 
10.00 





1970 10 1970 11 1970 12 


20.00 
20.00 
20.00 


6.67 
3.33 
10.00 


20.00 
20.00 


20.00 
20.00 
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2.3. THE USE OF OPERATORS IN LIST ARGUMENTS 


Many APL ECONOMETRIC PLANNING LANGUAGE functions such 
as TABULATE, PLOT, INTERSECT or REGRESS (see chapter 3) 
accept a LIST of rank one time series as input. It is 
permissible, where meaningful, to include operations in a 


LIST as demonstrated in the following example. 


Example 8 
2.4. EQUATIONS AND MODELS 


A single equation is made available to the system in 
the form of an APL-character vector containing a sequence of 
well-defined operations and an explicit result (dependent 
variable). As an example consider the first equation of the 
reference model: 


EQU1-'S1-(A11 T Wi D #2) P (A12 T 1 LAG 81) P 433 T DI! 
EQU1 ts the user-supplied name of the equation, Si the 
dependent variable. Intermediate results within the same 
equation are permitted, e.g. 


EQU1«'51-(A11 T WR-W1 D W2) PÅ 
A model [5 a set of equations. It [5 presented to the 
system in the form of a character matrix; each row contains 
one equation, A function 


NAME-DMODEL LIST 


assists the user in shaping this array. LIST contains the 
names of the equations. NAME is a user supplied name for the 
model. 


Example 9 


(Note that SALESMODEL is with the exception of a missing 
function name the canonical representation of a valid APL 
function. This fact will be utilized during the solution of 
the model, which in effect leads to a repeated execution of 
the corresponding APL function). 














+) 
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EXAMPLE 8 


1 1961 1 ST 0 


10 5 9 3 TABULATE 'S1,52,51 D 52! 


S1 


S1 D S2 


1961 1962 1963 1964 1965 


1613.578 1657.487 1731.530 1672.207 1702.119 
1971.696 1900.684 1808.007 1590.523 1633.202 
«818 «872 «958 1.051 1.042 
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EXAMPLE 9 @ 


EQU1-'S1-(227.2 T W1 D W2) 84) P .12 T DI' 


ኮህ 
ምጫ 
on 
ση 
InN 
= 
ኩ። 
I> 
Ic 


EQU2-'S2-( 80.5 T W1 D W2) 


ኮህ 


(.83 


m 

> 
Jt 
J> 
IQ 


S2) P .13 T DI' 


EQU3-'Wi1-.3 P (.34 T S1D 52) P ' 


EQU3-EQU3,' .15 T (1 LAG 51) D 1 LAG S2' 


EQUu-'W2-.39 P (7.29 T 


T 
EQUu<EQUu,',2 T (1 LAG S1) D 1 LAG 532! 


SALESMODEL1+DMODEL ‘EQU1,EQU2,FQU3,EQU4' 


SALESMODEL1 

S1<(227.2 T Wi D #2) P (.55 T 1 LAG 81) P .12 T DI 

S2<( 80.5 T #3 D W2) P (.83 T 1 LAG. 52) P .13 T DI 

W1+.3 P (.3u T 81 D S2) P .35 T (1 LAG S1) D 1 LAG S2 : 
W2+.39 P ( .29 T 81 D S2) P .2 T (1 LAG $1) D 1 LAG 82 | 








4) 
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2.5, THE USE OF OPERATOR SYMBOLS 


The APL ECONOMETRIC PLANNING LANGUAGE includes a 
translation capability, which permits the user to return to 
APL operator symbols for most of the mnemonic operator names 
and also to abbreviate LAG, DEL and RIO operators. For 
example a valid alternative for writing EQU1 in the 
preceeding paragraph would be: 


 EQU1-'S1-(A11xWi*W2)*(A12xS1c 12)*413xDI' 


However, this form cannot be directly processed by the 
functions for model solution (chapter h), but must be 
translated in advance. This task is performed by the 
function 


NAME2-TRANSLATE NAME1 


NAME1 15 the name of an equation or a model containing 
abbreviated operator symbols. WAME2 15 the name of the 
translation. 


Operations to be immediately executed can also be 
written in abbreviated form by means of the DO function 


( NAME+ ) DO A 


A ts a sequence of operations in character representation, 
e, g£. 


2-00 'X4Y' 
This ts equivalent to 
| 2+፲ P Y 
or 
DQ '2<X+Y' 


DQ flrst translates the right argument and then 
executes It. The last alternative formulation above 
demonstrates how the DO function can be used to solve an 


 equatlon (of a model), if the right hand side is completely 


specifled. But DQ does not accept a multi-equation model as 
Input. In this case the solution technique of chapter 4 must 
be applled. 


Table 3 gives a summary of all symbols permitted as 
abbreviations for mnemonic operator names. 
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TABLE 3: Valid operator symbols o 
mnemonic mnemonic 
operator symbo] operator symbol 
name name 

P + G > 

H x NE # 

T x AND A 

D + OR ν 

PW * NOT ~ 

L < K LAG X — Xc-K» 

LE = DEL X AX 

GE > RIO X OX 


As demonstrated in the following example, the ASSIGN 
function (see 1.1.) can also be used for the decomposition 
of a model Into its component equations. This capability [5 
advantageous if the equations of the final model form are to 
be stored separately on a file; see the section on DATA FILE 
ROUTINES, 


The last three abbreviations for LAG, DEL and RTQ can 
only be used if the right argument [5 a name of a serles but 
not a compound expression. Not permitted are, e.g., 

^ (X*Y) or (X*Y)c 1» 
These should rather be written: 

AX+AY or Xc 15+Yc 32 
or in long form: 


DEL X+Y or 1 LAG X+Y 
Example 10 
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E1-'91-(227.8xW1*W2)*(.55xS1c 12)*.12xDI' 


SALESMODEL2*-DMODEL 'E1,E£QU2,EQU3,EQU4' 


SALESMODEL2 _ 
S1+(227.4xW1:W2)+(.55xS1c 1>)+.12xDI 
S2<( 80.5 T Wi D W2) P (.83 T 1 LAG 52) P «19 T DI 
W1-.3 P (.34 T 91D S2) P .15 T (1 LAG S1) D 1 LAC 
W2+.39 P (7.29 T S1 D 82) P .2 T (1 LAG S1) D 1 LAG 
SALESMODEL3«TRANSLATE SALESMODEL2 
SALESMODEL3 
S1+(227.4 T W1 D W2) P (.55 T? (1 LAG S1)) P 
S2<( 80.5 T Wi D W2) P (.83 T 1 LAG S2) P «13 T DI 
Wi-.3 P (.34 T 51 D S2) P .15 I (1 LAG 51) D 1 LAG 
W2-.39 P (7.29 T S1 D S2) P .2 T (1 LAG S1) D 1 LAG 
'L1,L2,L3,Lut ASSIGN SALESMODEL3 
L3 
W1-.3 P (,35 T S1 D 52) P .15 T (1 LAG 81) D 1 LAG 


0 ST 1 1957 1 


DO 'Si+S2' 
495401 3446.646 3255.064 3358.905 3405.945 3410.322 


D 


Q 
195501 0. 


L3 | 
3357049 0.36697926 0.39319758 0.43409402 


S2 
S2 


22 
22 


22 


«12 


m 


DI 
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5. REGRESSION @ 


In section 3, several functions for the statistical 
estimation of the parameters of an econometric model are 
described. As mentioned earlier, it is assumed that the 
model is linear in its parameters. The following estimation 
techniques are provided, 


Ordinary Least Squares 
Generalized Least Squares, especially first order 
autoregressive disturbances 
Instrumental Variables Substitution (two-stage least 
Squares) 
Principal Components Substitution . 
Polynomial Distributed Lags (ALMON-Lags) 


The principal characteristics of these techniques will 
be briefly discussed in the subsequent chapters. For a 
rigorous introduction to the subject the reader is referred 
to econometric textbooks such as THEIL [8] or SCHNEEWEISS 
[7]. A concise treatment of the more technical aspects 15 
found in EISNER and PINDYCK [3] 





3.1. SIMPLE CORRELATION COEFFICIENTS (THEIL, p. 64) 


The pairwise simple correlation coefficients of a LIST 
of time series are calculated by 


( NAME- ) COR LIST 
Only data within the common intersection are included. Of 
course, this Intersection process can be explicitly 
controlled; see 1.5. Appropriate error messages are printed 
If the intersection process falls. It Should be menttoned 
again that COR and all other functions in this section 3 
work only for a LIST of rank one series, 


Example 11 


ti 


1 


=. 
0. 
Eo 
0. 


EXAMPLE 11 
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C-COR 'S1,92,W1,W2,DI' 


C 

70.9771 
9771 1 
9627 70.9531 
9121 0.9008 
8461 70.8972 





0.9627 

70.9531 
1 

0.9516 
0.7718 


0.9121 
0.9008 
~ 0.9516 
1 
0.6472 


0.8561 
0.8972 

0.7718 
~ 0.6472 
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3.2. ORDINARY LEAST SQUARES (THEIL, pp. 101) 


ለ1] parameter estimation techniques of the APL 
ECONOMETRIC PLANNING LANGUAGE must be initiated by calling 
the function REGRESS: ΕΘ 
ΝΑΜΕΕΙΥΙ REGRESS LIST 
WITH: 


The function responds with an inquiry WITH:, expecting the 
selection of the desired estimation technique. REGRESS 
continues to respond by means of the keyboard request WITH:. 
In that fashion several techniques can be combined within 
the same regression. A null carriage return finishes the 
selection phase. 


LIST is the list of independent variables (time 
series). The presence of the constant term is indicated by a 
"1", catenated to the list of variables; see example below. 
Y ts the dependent variable. The explicit result NAME [5 a 
character representation of the regression equation 
including the estimated values for the parameters. In this 
form the equation can be directly passed on to the 
construction of the complete model; see 2.4. 


REGRESS also prints several statistical measures. The 
series of disturbances can be obtained by calling the global 
variable U after the regression has concluded. If the 
disturbances are needed not only immediately but also at a 
later point a new name should be assigned to them, because 
the contents of U will be destroyed with the next following 
regression. The covariance matrix of the parameters is 


contained in the global variable B. 


Ordinary least squares estimation will be provided if 
the first WITH: ls answered by a null carriage return. Only 
those observations which lie in the intersection of all LIST 

series and the dependent variable Y are included in the 
regression. Again, of course, the intersection process can 
be controlled explicitly. This is demonstrated in the 
example where the first three years of the observation data 
will be excluded from the analysis . 


There are two functions to control] the output of 
REGRESS. For K<0 the command 


QUTPUT K 


suppresses the print-out of the statistics; for non-zero 
input X it provides a full print-out. For Χ-ο, the global 
variable B does not contain the covariance matrix of the 
parameters but an intermediate result. However the 
disturbances U are still computed correctly. 
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PRECISION K 


with a positive integer X<16, specifies the significant 
digits for the parameter values in the explicit result NAME 
of REGRESS, The QUIPUT and PRECISION specifications remain 
active also for subsequent regressions until they are 
respecified. A PRECISION SPEC ERR occurs for an invalid 


Input K, 


REGRESS itself might return the following error 
messages: REGRESS RANGE ERR |f the intersection of the 
observation data falls; REGRESS SING ERR if the intersection 
of the series contalned in LIST has less than full row rank, 
t.e. some of the series are linearly dependent. Further 
error messages may be incurred by the selected estimation 
techniques; see subsequent paragraphs. 


Example 12 
32.3. GENERALIZED LEAST SQUARES (THEIL, pp. 236) 


For ordinary least squares, it is assumed that the 
disturbances have a scalar covariance matrix in order to 
yleld consistent and efficient estimates: 


C-AxXI ልጸ scalar, J unit matrix. 


For generalized least squares techniques this assumption is 
relaxed: 


C-AxV A scalar, V symmetric positive definite matrix. 


The (rare) case, for which V is known in advance, can be 
handled by answering the WITH : inquiry with GLS V 


NAME-'Y' REGRESS LIST 
WITH:GLS V 
WITH: 


V is a symnetric, positive definite matrix with as many rows 
or columns as there are observations. An entry v(s,t) of V 
is, except for an unknown scalar factor, the covariance of 
the disturbances u(s) and u(t) of periods s and t. 


A GLS DIM ERR ls returned if the dimension of V is not 
properly chosen. GLS PD ERR indicates that V is not 
symmetric and positive definite. 
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PRECISION 6 


1 1956 1 ST O 


5 ሀሀ4+፻2፤ REGRESS '1,S1 D S2,(1 LAG S1) D 1 LAG $2! 


WITH: 


COEF/VALUE/ST ERR/T-STAT..... 


1 .38506 ~O4656 8.26945 
2 .29101 .19820 71.46827 
3 . 19868 .16595 1.19725 


NO OF VARIABLES......... 2.00000 
NO OF OBSERVATIONS...... 10.00000 
SS DUE TO REGRESSION.... «00204 
SS DUE TO RESIDUALS..... «00352 
F-STATISTIC............. 1.99692 
STANDARD ERROR.......... . 02213 
R*2 -STATISTIC.......... . 36328 
R*2 CORRECTED........... . 18136 
DURBIN WATSON STATISTIC. 1.13518 


ጀፀሆዛ 
ሆ2ሩ ( 0.385065 7 1 


) 0.291008 T S1 D S2 ) P ( 0.198681 
T (1 LAG S1) D 


ጅ ( 
4 LAG 52) 
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e | a FOR THE SAKE OF SIMPLICITY IT MAY BE SOMETIMES 
a CONVENIENT (BUT NOT NECESSARY) TO DEFINE 
a NEW SERIES FOR THE COMPOUND EXPRESSIONS: 
PAS | | 
SR-S1 D 22 

: SLi«1 LAG S1 
SL2+1 LAG S2 
SRL*SL1 D SL? 


WR+W1 D W2 


1 1956 1 ST O 


EQU1+'S1' REGRESS 'WR,SL1,DI' 
WITH: 





COEF/VALUE/ST ERR/T-STAT..... 


1 227.17691 78.43286 2.89645 


2 54654 .06132 8.91224 

3 .12025 .03609 3.33226 
NO OF VARIABLES......... | 8.00000 
NO OF OBSERVATIONS...... 10.00000 
SS DUE TO REGRESSION.... 23254099.92031 
SS DUE TO RESIDUALS..... 13377.08146 
Ρ-6ΤΑΤΙΘΤΙΟ.........»... 1056, 15878 
STANDARD EFRROR.......... 43.71512 
R*2 -STATISTIC.......... .99943 
R*2 CORRECTED........... .99918 
DURBIN WATSON STATISTIC. | 2.05822 

EQU1 


Sie ( 227.177 T WR ፻ ( 0.546535 T SL1) P ( 0.120248 T DI ) 
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More important and more frequently used Is the 
specific case of a non-scalar covartance matrix, where the 
disturbances are generated by a first order autoregressive 
process (THEIL, pp. 250) 


U-S P RHO T 1 LAG U 
with 1<RHO<1 a scalar constant and Š a random sertes with 
zero expectation and scalar covariance matrix. This case [5 
treated by means of the command 41 RHO or 41 10 if RHO is 
not known [ከ advance: 


NAME-'Y' REGRESS LIST 
WITH:A1 RHO (10) 
(RHO:NN) | 
WITH: 
[የ RHO=10, an estimate for RHO [5 computed (THEIL, p. 254) 
and printed (WN) before the next WITH response, 


An 41 VAL ERR occurs if 1s|REO, A REGRESS SING ERR may 
be detected before RHO ls calculated; see 3.2. 


It should be noted that after the application of 
generalized least squares the transformed disturbances P*.xU 
are used in the evaluation of the statistics. P results from 
a factorization of the inverse of V: 


((QP)+.xP)=8V 


In the case of the first order autoregressive process this 
corresponds to the use of the disturbances Š Instead of the 
original U, Therefore, after the regression the global 
variable Ü contains the series $5, | 


Example 15 


3.4, INSTRUMENTAL VARIABLES SUBSTITUTION (two-stage least 
squares) | 
(THEIL, p445, pp. 451) 


One of the basic assumptions for ordinary least 
Squares to yleld unblased, consistent and efficient 
estimates Is, that the right hand side variables in LIST are 
uncorrelated with the disturbance term. This assumption | 
cannot be maintained when the regression equations are part 
of a simultaneous equations model. For the reference mode] 
it can easlly be seen that tn the first equation W1 D ሆ2 
must be correlated with U1, stnce equations 3 and 4 contain 
on. the right side S1 and therefore also the disturbance term 
U1. | 
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 &X4MELE 13 


0 ST O 
OUTPUT 0 


PRECISION 4 


a GENERATION OF A (12,12) RANDOM COVARIANCE MATRIX 


a WHICH WILL BE USED FOR THE GLS DEMONSTRATION: 


VTEST+12 12pRECT 12x12 


VTEST+VTEST+.x&VTEST 





GEQU3+'W1' REGRESS '1,SR,SRL' 
WITH: GLS VTEST 
WITH: 


GEQU3 ` _ 
Wie ( 0.01531 T 1 ) P ( 1.523 T SR ) P ( 70.9761 T SRL) 


QUTPUT 1 


PRECISION 6 


&1 


J AEQU3-'W1' REGRESS '1,SR,SRL' 
WITH: Δι 10 
RHO: 0,2297705705. 
WITH: 


1 «07880 «05206 1.51348 


2 1.01861 .30507 3.33890 
3 .60711 «28230 2.15055 


NO ΟΕ VARIABLES......... 3.00000 
NO OF OBSERVATIONS...... 12.00000 
SS DUE TO REGRESSION.... 1.43788 
SS DUE TO RESIDUALS..... . 01793 
P-STATISTIC............. 240.62091 
STANDARD ERROR.......... . 04463 
R&2' -ΦΤΑΤΙΟΤΙ͂Ο....».».... .98769 
R*2 CORRECTED........... . 98358 
DURBIN WATSON STATISTIC. 1.65762 
AEQU3 


Wi« ( 0.0787991 T 1 ) P ( 1.01861 T SR ) P ( 70.607108 T SRL 
) . 
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In ordec to remove this deficiency one can apply the 
technique of instrumental variables substitution. The basic 
idea behlnd this technique is to regress all variables in 
LIST which are assumed to be correlated with the disturbance 


term, first on a set of auxlllary variables (Instrumental 


varlables), before the malin regression [5 performed. These 
auxl11ary varlables should on the one hand be strongly 
correlated with the original variables, on the other hand ` 
uncorrelated with the disturbance term. The results of these 
"first stage" regresslons replace the original variables in. 
LIST, Since the substitutes are linear combinations of the 
Instruments they are also uncorrelated with the disturbance 
term. The main regression can now take place with the 


modtfted LIST contents. 


It had been shown independently by BASMANN and THEIL, 
that such a two stage procedure ylelds asymptotically ` 
unblased, consistent and efficient estimates if all 
exogenous serles of the complete model are used as a set of 
Instruments (THEIL, pp. 497). In that case the procedure had 
been given the name two-stage least squares technique. 


| Often, at least for very. large models, the set of 
exogenous variables [5 quite large and some serles are 
htghly collinear, so that a full two-stage least squares 
procedure would either fall or at best yleld inaccurate 
results. A popular way out Is the use of only a subset of 


the Instruments. Here, however, a problem of selecting the 


best subset arises. This is partly resolved by the use of 
principal components of the total set. The latter technique 
will be discussed in the next paragraph. 


Instrumental variables substitution Is effected by 
means of the command INST 


NAME-'Y' REGRESS LIST 
WITH: ILIST INST LISTA 
WITH: . 


ILIST ls the list of instruments, LISTA a sublist of LIST, 
contalning the series for which the substitution is to be 
performed. Further INST commands can be added if different 
subsets from LIST should be regressed on different 
instrumental variables: 


WITH: ILIST1 INST LISTA 
WITH:ILIST2 INST LISTB 
WITH: 
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The instruments must have the same periodicity and the 
same number of observations as the original series in LIST 
(after intersection). Otherwise an INST RANGE ERR wil] 
occur. An INST SING ERR ls delivered if some of the 
instruments are linearly dependent. An INST SPEC ERR 
indicates that the right argument of INST ts not a subset of 
LIST. 


The following example demonstrates a full two-stage 
least squares estimation of equation 2, where all exogenous 
variables will be used as Instruments. | 


Example 14 


3.5. PRINCIPAL COMPONENTS SUSTITUTION (THEIL, pp. 4&6) 


This technique provides the user with a systematic 
procedure of selecting the right Instruments If the complete 
set of exogenous variables Is too large. The basic tdea ts 
to transform the complete, normaltzed set of Instruments 
into a set of as many new series, principal components, 
which have two properttes: They are orthogonal, thus 
avoiding the problem of collinearity, and they are ordered | 
In a sequence according to thelr decreastng contribution to 
the explanation of the original series (in terms of variance 
decomposition). Finally, as many principal components are . 
selected to be instruments as are needed to exceed a LE 
user-specified minimum variance level. The specified . REUS 
variance level should be less than the total variance, 
because the selection of all principal components is 
equivalent to the immediate use of all exogenous variables 
as instruments and no advantage will be galned. 





Principal components substitution can be performed by x 
the PCINST command | | 


NAME-'Y' REGRESS LIST 
 WITH:ILIST PCINST LISTA 
VARIANCE:K 
WITH: 


As in the previous chapter, ILIST ls a set of instruments. 
Here they will, as a first step, be transformed Into 
principal components (ustng the JACOBI method for 
calculating eigenvectors). LISTA 15 a subset of LIST, for 
which principal components substitution [5 to take place. 
PCINST responds with VARIANCE:. The answer X [5 elther a 
positive 
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A EXAMPLE 14 


IEQU2-'S2' REGRESS 'WR,SL2,DI' 
WITH: 'SL1,98L2,1,SRL,DI' INST 'WR' 
WITH: | 


COEF/VALUE/ST ERR/T-STAT..... 





1 32.69828 103.63321 343558 
2 «86337 .05982 14.43311 
3 .04115 .08363 49205 
NO OF VARIABLES......... 3.00000 
NO OF OBSERVATIONS...... 12.00000 
SS DUE TO REGRESSION.... 53964323.76781 
SS DUE TO RESIDUALS..... 74063.20131 
P-STATISTIC............. 2185.87596 
STANDARD ERROR.......... 90.71518 
R*2 -STATISTIC.......... .99863 
R*2 CORRECTED........... .99817 
DURBIN WATSON STATISTIC. 1.72067 
IEQU2 - 
82+ ( 32.6943 T WR ) P ( 0.863371 T SL2) P ( 0.0411482 T DI 
) 


ILIST-'1,SL1,SL2,SRL,DI' 
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real number less than 1, or integral. In the first case, K 
denotes the percentage of total variance to be explained. In 
the second case, K is the number of selected principal 
components, the variance criterion being ignored. This 
number should be smaller than the number of series contained 
in LIST, not counting a "1", Again several principal 
components substitutions can be performed within the same 
regression; see 3,4, | 


A PCINST RANGE ERR is returned if the instruments in 
ILIST do not have the same periodicity and the same number 
of observations as the series in LIST (after the 
intersection process). A PCINST SING ERR occurs if some of 
the instruments in ILIST are linearly dependent. A PCINST 
SPEC ERR ts delivered if LISTA ts not a subset of LIST, 


Example 15 
5.6. POLYNOMIAL DISTRIBUTED LAGS (ALMON [11 ) 


An equation may incorporate a welghted sum of lags of 
a series X | | 


(40 T X) P (43 Z ፲ 1>) P ..... (AK T Xc-K»2) 


A direct attempt to estimate this term, possibly in 
combination with other terms in the same equation, may fall 
because of collinearity or insufficient degrees of freedom. 
A common approach for resolving this deficiency ts to assume 
a-priori that K+1 parameters are determined by a polynomial 


of degree less than K, It is also common to impose constraints 
on the shape of the polynomial, as desired. This usually takes 


the form of constraining the beginning pun and/or end (tail) 
of the polynomial to be zero (see below). | 


Polynomial distributed lags are performed by a DIST 


command following a WITH response. The syntax is 


NAME-'Y' REGRESS LIST 
WITH:N DIST 'X' 
WITH: 


X is the member of LIST for which distributed lags are to be 
estimated. LIST contains only the unlagged term; the correct 
lags are provided automatically at the end; see example 

below. Within each equation one or more series can be expanded 
into a weighted sum of lags. 
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A EXAMPLE 15 


PCEQUu-'W2' REGRESS '1,SR,SRL' 
WITH: ILIST PCINST 'SR' 
VARIANCE: .95 
WITH: 


1 «50882 .05640 9.02110 
2 7.94805  .38712 2.54897 
3  .75369  .35350 2.13208 


NO OF VARIABLES......... 2.00000 
NO OP OBSERVATIONS...... 12.00000 





SS DUE TO REGRESSION.... .01501 

SS DUE TO RESIDUALS..... ,00928 

FPSOPATISTIC 3406606508054 7,27822 

STANDARD ERROR.......... «03211 

R*2 -STATISTIC.......... .61794 

፳#2 CORRECTED........... .53304 

DURBIN WATSON STATISTIC. .95153 
PCEQU4 

W2+ (0.508825 7 1) P ( 0.948046 T SR ) P ( 0.753688 T SRL 
) 





equal to the number of variables. Note, that the effective 
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The BIZZ command can be called repeatedly during the same re- 
gression for different independent variables. o 


N is a three-element vector. ΔΙ] supplies the 
number of desired terms in the lag structure for x (i.e., 
K+1), N[2] is the degree of the polynomial, and #[31] in- 
dicates the constraints imposed: 0 = no constraints, 1 = 
zero-constraint on head of polynomial, 2 = zero-constraint 
on tail, 3 = zero-constraint on both head and tail. For 
the special (common) case where N[2 3] = 3 3,these two 
elements may be omitted and only N[1] need be typed, 


The original series,X is used for the construction of 
the lags, regardless of a previous generalized least squares 
or an instrumental variables transformation; see 3.7. If X 
has not enough values in front of the current intersection, 

a reduction of this intersection may occur. A DIST RANGE ERR 
is produced if the number of periods becomes less than or 


number of variables is determined by the original contents 
of LIST and the specified degree of the polynomial, taking 
the specified constraints into account. A DIST SPEC ERR in- 
dicates a misspelling of one of the arguments of DIST. 


Example 16 
3.7. COMBINATION OF ESTIMATION TECHNIQUES e 


From a purely technical point of view the user may 
combine within the same regression different estimation 
methods in an optional fashion. In the example below, the 
correction of the disturbance term and an instrumental 
variables substitution are applied to the same equation. 
All techniques are executed in the same sequence as they 
were selected. 


However, the chosen sequence affects the statistical 
properties of the estimates. The theory of correct selection 
is not yet fully explored. Some aspects are discussed in 
EISNER and PINDYCK [3]. They suggest the following sequence: 


Generalized least squares (GLS, A1), 


instrumental variables substitution (INST, PCINST), 


polynomial distributed lags (DIST), 


ordinary least squares. 


Example 17 o 
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DLEQU2-'S2' REGRESS 'WR,SL2,DI' 
| WITH: 4 2 0 DIST 'pT! 
WITH: 





1 ~227.48477 216.41694 — “1.051185 
2 0.80948 0.07789 10.39235 
3 "0.34414 0.53270 70.64602 
4 70.29938 0.42354 70.70684 
5 0.08844 0.47402 0.18657 
6 0.81930 0.46762 1.75207 
ΝΟ OF VAPIABLES......... 5.00000 
NO OF OBSERVATIONS...... 10.00000 
SS DUE TO REGRESSION.... 38721690.46784 
SS DUE TO RESIDUALS..... 35019.72157 
PeSTATISTIC.ISae ολο መ 1105.71098 
STANDARD ERROR.......... 83.68957 
፳#ቁ2 -STATISTIC.......... 0.99910 
R*2 σΟΒΕΕΟΤΕΡ........... 0.99819 
DURBIN WATSON STATISTIC. 2.64357 
 DLEQU? 
S24 (7227.485 T WR ) P ( 0.809484 T 652 ) P ( 0.344138 T DI 
) P ( 70.299375 T 1 LAG DI) P ( 0.0884364 T 2 LAG DI) P ( 0 


«819297 T 3 LAG DI) 


& 9 


1 1956 1 ST O 


MEQU2+'S2' REGRESS 'WR,SL2,DI! 
WITH: A1 10 
RHO: 0.1972740522 
WITH: ILIST INST 'WR' 
WITH: 


COEF/VALUE/ST ERR/T-STAT..... 


1  87.07280 162.74288 53503 

2 .82908 .07589 10.92410 

3 .12147 .09732 1.24806 
NO OF VARIABLES......... 3.00000 
NO OF OBSERVATIONS. ..... 10.00000 
SS DUE TO REGRESSION.... 27630848.94413 
SS DUE TO RESIDUALS..... 63664.82282 
F-STATISTIU. ¿us ውሬ 1012.67824 
STANDARD ERROR.......... 95.36758 
R*2 -STATISTIC.......... .99770 
Rx2 CORRECTED........... .99672 
DURBIN WATSON STATISTIC. 1.65174 

MEQU2 


S2< ( 87.0728 T WR ) P ( 0.829077 T 652) P ( 0.121566 T DI 
) 
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MEQU3-'Wi1' REGRESS '1,SR,SRL' 
WITH: A1 10 
RHO: 0.1732689385 
WITH: ILIST INST 'SR' 
WITH: 


1 . 30444 .05852 5.20197 


2 .33337 «26026 1.28091 
3 .15827 «21698 -, 68333 


NO OF VARIABLES......... 3.00000 
NO OP OBSERVATIONS...... 10.00000 
SS DUE TO REGRESSION.... 2.98019 
SS DUE TO RESIDUALS..... . 00508 
F-STATISTIC............. 1368.46052 
STANDARD ERROR.......... . 02694 
R*2 -W-STATISTIC.......... . 99830 
R*2 CORRECTED........... «99757 


DURBIN WATSON STATISTIC. 1.98565 
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MODEL SOLUTION 
SOLUTION OF A SIMULTANEOUS EQUATIONS MODEL 


A simultaneous equations model which may be nonlinear 
in its variables is solved in the APL ECONOMETRIC PLANNING 
LANGUAGE by means of the GAUSS-SEIDEL-technique; see 
ORTEGA and RHEINBOLDT [6]. The essential idea of this 
iterative technique is to first normalize the equations 
of the model in such a way that each unlagged endogenous 
variable (solution variable) appears on the left hand side 
of one and only one equation. An initial guess for the 
first period solution of the left hand side variables being 
supplied, the model equations are then successively evaluated 
in order to yield new and improved estimates for the solu- 
tion. Once a solution for a period is found, the algorithm 
takes this solution as initial guess for the next period, 
and the iteration process starts again for that period. 

The k-th iteration would be (period index omitted): 


Yl(k)« Fl (Y2(k-1),Y3(k-1),...,Yh (k-1)) 
Y2(k)+ F2(Y1(k), Y3(k-1),...,Yn (k-1)) 


YN(k)« Fn (Yl(k), Y2(k) ,...,Yn-1(k)) 


The right hand side of equation n may contain the series 
Yn only in lagged form, i.e. the term K LAG Yn with Και: 


The criterion for convergence is 


EPS G ABS (Xn(k) M Yn(k-1)) D Yn(k-1) 


for all series Yn. EPS is a predefined, sufficiently small 
positive scalar. 





dr; 
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@ If the convergence process is oscillatory, a relaxation parameter 
ALPHA might provide damping: 


| Yn(k)<(ALPHA T Fn(k))P (1-ALPHA) T Yn(k-1) 
0<ALPHAS1 


| where Fn(k) is the right hand side of the unrelaxed process above. 
| The solution process is initiated by the function SOLVE: 
LIST SOLVE NAME 


| NAME is the name of the model. As mentioned earlier the model 
| must be presented in normalized form. This is not done automatically 
but has to be considered by the user when defining his equations. 

It is good practice to specify the model from the very beginning 

in normalized form, which is always possible by the addition of 
appropriate definition equations. LIST is the list of all exogenous 





series in the model, excluding lagged endogenous series. 


Before SOLVE can be requested, several preparatory steps are 
necessary. First, the simulation horizon, i.e. the periods for 
which the model should be solved, is to be specified via the setting 
of a time frame; see 1.5. The horizon may extend over past periods 
for model validation (ex post simulation) or over future periods 
for the purpose of forecasting (ex ante simulation). 


| Secondly, all exogenous series must be available within the 
| active workspace (see e.g. the file routines in chapter 5 as one 
| example to load series). The series must have the proper length 
| in order to fill the model requirements for yielding solutions over 
| the whole simulation horizon; variables appearing in lagged form 
| must include pre-horizon values depending on the lag structure. 
| Also for lagged endogenous series pre-horizon values in proper 
e length must be present. 
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Initial guesses for the endogenous series within the simulation 


horizon are computed automatically. Only for the first period, 

a predefined value is observed. For subsequent periods the solutions 
of the period before are taken. If no initial value for the first 
Simulation period of an endogenous series is found, the value from 
the closest period to that is taken (either the last or the first 
value of that series depending of its time structure). If the 
series is not defined at all, the square root of 2 is assumed as 
initial value. 


Finally four different parameters are to be set to control the 
solution process and output 


SOLVECON EPS 


sets the convergence bound EPS. A typical value is EPS<.001 
(O.1 percent relative error). 


SOLVEIT K 


defines a limit K for the number of iterations per period. 
SOLVEREL ALPHA 


sets the relaxation parameter ALPHA. If it is assumed that the 
solution process does not oScillate from iteration to iteration 
then ALPHA<1 is appropriate, otherwise a value between zero and one 
should be specified. 


Output is controlled by means of 
QUIPUT K 
(compare chapter 3.2). For K«0 only a convergence message per period 


including the number of required iterations is printed, otherwise 
a full convergence diagnostic is produced; see example 18. 
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The control parameters remain active for repetitive applications 
of SOLVE until they are respecified explicitly. They can also be 


key, execution of SOLVECON, SQLVEIT, SOLVEREL or QUTPUT, and return 


to the interrupted line of SOLVE. 


respecified during execution of SOLVE by interruption with the ATTN 


After completion of SQLVE the solution series are available 
to the user by the names given in the model; see example 18. If 
these solutions are to be compared with the observation data for 
the same variable, their names should be respecified to distinguish 
between equally named series in the model and on the data file. 


If convergence is not achieved within the permitted iteration count, 


SOLVE terminates with a message ITERATIOK LIMIT EXCEEDED. The 


endogenous series then contain the "solutions" found so far. 


A message MISSING LAGS IN EQU Yn is delivered if not enough 
pre-horizon values for a lagged series in that equation are supplied. 
An EQU N ERR occurs if the n'th equation contains a syntax error 
which prevents the APL execution of the model (internally the model 
is transformed into an APL function with the name CMR). An EXOG 
RANGE ERR indicates that one or several of the exogenous series are not 
properly specified. An APL length error during execution of CMR may 
result from the omission of an exogenous series in the left argument 
LIST of SOLVE. Various intersection or model language errors occur 


if errors in the model equations are detected during execution. 


Example 18 
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8 ASSEMBLY OF THE MODEL FROM THE REGRESSION RESULTS 
ϐ (SECTION 3), ADDING TWO NEW EQUATIONS FOR P1 AND P2 


EQU1 | | | 
S1<(227.177 T W1 D W2) P (.546535 T 1 LÆ δι) P (.1202u8 T DI) 


MEQU2 


S2<( 87.0728 T Wi D W2) P (.829077 T 1 LAG S2) P (.121466 T DI 


) 

MEQU3 B 
W1i«(.304443 T 1) P (.333374 T S1 D 82) P (^ .158266 T (1 LAG 8S1 
) D 1 LAG S2) 


| EQUN | | 
W2+ ( 0.385065 Z 1 ) P ( 0.291008 T 81 D 82 ) P ( 0.198681 T 
(1 LAG 51) D 1 LAG 52) | 


EQU5<'!P1<S1 D W1' 
EQU6+'P2+S2 D W2' 


EQULIST+'EQU5,EQU1,MEQU3 ,EQU6 ,MEQU2 ,EQUY' 


SALESMODEL*DMODEL EQULIST 


SALESMODEL 
P1<S1 D Wi 


S1+(227.177 T Wi D W2) P (.546535 T 1 LAG #1) P (.120258 T DI) 


ሆ4+(, 305553 T 1) P (,333375 T S1 D S2) P (7.158266 T (1 LAG S1 


) D 1 LAG S2) 
P2-82 D W2 


= 
it 
m 
ΚΩ 
ta 
N 
ሠ 
ህዝ 


S2<( 87.0728 T W1 D #2) P (.829077 T (.121466 T DI 
) | | 


T 
P ( 0.291008 T S1 D S2 ) 


W2+ ( 0.385065 T 1 ) P ( 0.198681 7 
(1 LAG S1) D1 LAG S2) | | 
ο ST 0 
P1+DO EQUS 
P2-DO EQUO6 
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A SIMULATION HORIZON 


11954 1 ST 0 


a I. F. THE SIMULATION COMPRISES EX POST THE WHOLE 
a OBSERVATION PERIOD EXCLUDING THE FIRST YEAR 1953, 


@ WHICH IS NEEDED TO FILL THE LAGS IN THE MODEL 


A SPECIFICATION OF CONTROL PARAMETERS 


SQLVECON .001 
SOLVEIT 100 
SOLVEREL 1 
QUTPUT 1 
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ο EXECUTION OF SOLUTION ALGORITHM 


'DI' SOLVE SALESMODEL 
PER 1 IT 1 MAX DEV 1.805301706 AT VAR Wi 
NON CONV VARS: 
S1,/1,92,W2 
PER 1 IT 2 MAX DEV 0.5354367893 AT VAR P1 
NON CONV VARS: 
P1,91,/91,P2,S92,W2 
PER 1 IT 3 MAX DEV 0.2504780124 AT VAR P1 
NON CONV VARS: 
P1,51,W1,P2,S2 
ATTENTION 
SOLVEL48] 
QUTPUT 0 
>48 
PER 1 1954 1 CONVERGED AT ITERATION 6 


PER 1 1955 1 CONVER ED AT ITERATION 6 
PER 1 1956 1 CONVERGED AT ITERATION 6 
PER 1 1957 1 CONVERGED AT ITERATION 6 
PER 1 1958 1 CONVE® ED AT ITERATION 6 
PER 1 1959 1 CONVERGED AT ITERATION 6 
PER 1 1960 1 CONVERGED AT ITFRATION 6 
PER 1 1961 1 CONVERGED AT ITERATION 7 
PER 1 1962 1 CONVERGED AT ITERATION 7 
PER 1 1963 1 CONVERGED AT ITERATION 7 
PER 1 1964 1 CONVERGED AT ITERATION 7 


PER 1 1965 1 CONVERGED AT ITERATION 4 
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በ SOLUTIONS 


51 
195501 711.9592266 1002.738006 1191.858813 1312.992263 
1393.424057 1469.837325 1533.138117 1593.946902 1658.992984 
1720.965193 1789.765572 1853.13183 


S 2 
195501 2869.002355 2610.158761 2402.118477 2226.224512 
2073.280515 1953.238 1851.484435 1767.282634 1702.933991 
1652.393083 1626.219559 1625.935757 


ES1+S1 
ፐ22+22 
EW1<W1 
HW2<W2 





SAVED 


S1 


ES1 


22 


ES2 


W1 


EW1 


W2 


EW2 





"COMPARISON 


1954 


398.658 


711.959 


2856.406 


2869.002 


0.126 


0.376 


0.425 


0.328 


WITH 


5 6 9 3 TARULATE LIST 


1955 


690.081 


1002.738 


2668.824 


2610.159 


0.236 


0.396 


0.371 


0.323 
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)COPY EPLDEMO BASSDATA 
15.53.13 03/11/76 


1956 


952.423 


1191.859 


2453.522 


2402.118 


0. 367 


0.413 


0.360 


0.317 


THE OBSERVED DATA 


1957 


1232.111 


1312.992 


2178.211 


2226.225 


0.453 


0.427 


0.279 


0.312 


LIST*'S1,FS51,,52,FS2,,W1,FW1,,W2,EW2,,,' 


1958 


1414.718 


1393.424 


2008.075 


0.451 


0.441 


0.289 


0.307 


1959 
1526.195 
1469.837 


2006.827 


1953.238 


0.495 


0.456 


0.287 


0.300 











22 


ES2 


W1 


EW1 


W2 


EW 2 


1960 
1571.303 
1533.138 


1994930 


1851.484 


0.462 


0.469 


0.294 


1961 
1613.578 
1593.947 


1971.696 


1767.283 


0.455 


0.482 


0. 286 


0.287 
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1962 


1657.487 


1658.993 


1900.684 


1702.934 


0.434 


0.495 


0.283 


0.281 


1963 
1731.530 
1720.965 


1808.007 


1652.393 


0.518 


0.507 


0.285 


0.276 


1964 


1672.207 


1789.766 


1990. 523 


1626. 220 


0.506 


0.517 


0.283 


0.272 


1965 


1702.119 


1853.132 


1633.202 


1625.936 


0.486 


0.521 


0.317 


0.272 
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35 50 10 PLOT 'S1,ES1' 


2000- 
1750- 
1500- 
1250- 
1000- 

750- 

500- 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
50-| 
| 
0 


0. 245 5.0 Fes 10.0 
ABSCISSA = TIME STARTING FROM 1954 
ο = $1 


* 


ES1 


12,5 











PER 1 


PER 1 


PER 1 


S1 


22 


W1 


W2 


DI 


59.1 


e PRODUCING A FORECAST 1966 TO 1968 
e FORECAST OF EXOGENOUS VARIABLE DI: 
DI<DI,3800 4050 4250 

n SIMULATION HORIZON 

1 1966 1 ST 1 1968 1 

n CONTROL PARAMETERS REMAIN UNCHANGED 
'DI' SOLVE SALESMODREL 

1966 1 CONVERGED AT ITERATION 10 

1967 1 CONVERGED AT ITERATION 8 


1968 1 CONVERGED AT ITERATION 7 


40 3 9 3 TABULATE 'S1,S2,,W1,W2,,DI' 


1966 1967 1968 


1821.636 1902.050 1951.880 


1649.042 1698.406 1770.613 


0.518 0.514 0.506 


0.271 0.279 0.287 


3800.000 4050.000 4250.000 


60 


4.2. REORDERING THE SEQUENCE OF THE EQUATIONS OF A MODEL መ 


The APL ECONOMETRIC PLANNING LANGUAGE contains a 
function ORDER which applies VAN DER GIESSEN‘s algorithm 
[9] to reorder the sequence of equattons of a model. Its 
aim is an almost recursive model, which contains on the 
right hand side as few variables as possible that were not 
already specified on the left hand side of previous 
equations. Such a model representation has the advantage 
that the most recent estimates can be used for each | 
iteration during the solution process. Thus, convergence can 
be accelerated. The reordering algorithm, furthermore, finds: 
a maximum list of endogenous variables which need not be 
assigned initial estimates because those are automatically 
calculated during the first iteration. 


The reordering is performed by 


NAME2+LIST ORDER NAME1 

NAME1 is the name of the unordered, NAME2 the name of the 

ordered model. By means of the LIST argument the user can 

select equations to be excluded from the reordering process. 

They are placed In the front of the reordered model, LIST 

contains the names of the endogenous series, f.e., the 

equations with those names as left hand side varlables. e 
LIST-' ' tf all equatlons should participate In the 

reordering algorithm. 


In addition to the ordered model (as explicit result), 
ORDER returns two lines of information: The sequence of the 
equations in terms of the left hand side variables and the 
list of variables for which initial estimates have to be 
provided hefore the solution process is entered. 


ORDER returns an ORDER SPEC ERR message 1f the names 
in LIST cannot be identified as left hand side variables. 








60.1 


A function 





| STRUCTURE NAME 


| displays the causal structure of a model NAME in graphic form. 
| The appearance of endogenous variables in the various equations 
| is marked by asterisks. The number of variables with entries 

| above the diagonal is an indication for the degree of recursive- 
| ness of the model in the assumed order of equations. 


| Example 19.1 
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SALESMODEL 
P1+S1 D W1 


S1<(227.177 T Wi D W2) P (.546535 T 1 LAG S1) P (.120248 T DI) 


1) P (.33337u T 81 D S2) P ( .188266 T (1 LAG δι 


W1+(.304443 T 


) D1 LAG 52) 
P2<S2 D W2 


S2<( 87.0728 T W1 D #2) P (.829077 T 1 LAG S2) P (.121466 T DI 


) 
/2ሩ ( 0.385065 T P ( 70.291008 T S1 D 82 ) P ( 0.198681 T 


T 1 ) P 
(1 LAG S1) D 1 LAG S2) | ; 


SALESMODORD<' ' ORDER SALESMODFL 
SEQUENCE OF EQUS: 
015,024572141,;P2,W2 


GIVE INIT APPR FOR: 
W1,W2 


SALESMODORD | 
$1-(227.177 T W1 D W2) P (,5ዛ6535 T 1 LAG 21) P (.120248 T DI) 


82-( 87.0728 T W1 D W2) P (.829077 T 1 LAG S2) P (.121466 T DI 


) 
P1<S1 D W1 


P, 


Wi«(.304443 T 1) P (.33337 T S1 D 82) P (7.148266 T (1 LAG S1 


) D 1 LAG 52) 


P2<S2 D W2 i 
W2+ ( 0.385065 T 1) P ( 70.291008 T 81 D 82 ) P ( 0.198681 T 
(1 LAG S1) D 1 LAG S2) 
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TRUCTURE SALESMODEL 


EQUATIONS 
ENDOG VARIABLES 


S 
P SW 
1 1 1 
P1| ο * x 
S1 | ο x 
W1 | * ο 
P2 | 
S2 | * 
W2 | * 
PSW 
1 1 1 
ROWS = 
COLUMNS = 
STR 
S SP 
121 
S1| ο 
S2 | o 
P1 κ ο 
Wil κ κ 
P2 | * 
W2| * * 
ο. 
1 2 1 
ROWS = 
COLUMNS = 


EQUATIONS 
ENDOG VARIABLES 


SAVED 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 


PER 1 
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)COPY EPLDEMO BASSDATA 


15.53.13 03/11/76 
ο ST 0 
P1<DQ EQU5 
P2+D0 EQU6 

1 1954 1 ST 1 1965 1 

'DI' SOLVE SALESMODORD 
1954 1 CONVERGED AT ITERATION 
1955 1 CONVERGED AT ITERATION 
1956 1 CONVERGED AT ITERATION 
1957 1 CONVERGED AT ITERATION 
1958 1 CONVERGED AT ITERATION 
1959 1 CONVERGED AT ITERATION 
1960 1 CONVERGED AT ITERATION 
1961 1 CONVERGED AT ITERATION 
1962 1 CONVERGED AT ITERATION 
1963 1 CONVERGED AT ITERATION 
1964 1 CONVERGED AT ITERATION 
1965 1 CONVERGED AT ITERATION 
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61.3 


SINGLE EQUATION SOLUTION 


For the ex post analysis and validation of a simultaneous 
equations model it may be desirable to solve each equation 
separately and independently from the results of other equations. 
Thereby, all right hand side variables are considered exogenous 
and must be completely specified before execution. In general, 
observation data are taken for all these variables rather than 
model solutions for part of them, as it would be in the simul- 


taneous equations solution. The single equation solution process 


is straight forward and consists, with the exception of feeding 
back the dependent variable (see below), of one execution per ` 
equation. This is comparable to the application of DO(see 2.5). 


‘Phe function 
SSOLVE NAME 
performs single equation solution of a model with a specified 


NAME without feedback. Feedback in that context means, that if 
the left hand side dependent variable is contained in lagged 


form also on the right hand side of the same equation, the 


equation result is substituted for the lagged series. Feedback 
can be applied by calling 


SSOLVEF NAME 
instead of SSOLVE. NAME is again the model name. 


After execution of either of the two single equation 
routines the solutions are available by the names given in the 


model. With SSOLVE some of the returned solutions may be 


shorter because of diadic operations on lagged series in the 
respective equations. | 


Example 19.2 
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)COPY EPLDEMO BASSDATA 
SAVED 9.31.06 01/22/76 


EQU1<!S1<(227.177 T W1 D W2) P (.546535 T 1 LAG S1) P ' 
EQU1+EQU1, '(.120248 T DI)! 
EQU2-'S2«(7 87.0728 T W1 D W2) P (.829077 T 1 LAG S2) P ? 
EQU2+EQU2,'(.121466 T DI)! 


EQU3<'W1<(. 304443 T gt P (.333374 T S1 D S2) P (7.148266 ! 
EQU3-EQU3,'T 1 LAG S1 D S2)' 


ኛ9ሆ4+፤ W2«(. 385065 T 1) P (7.291008 T S1 D S2) P (.198681 ' 
EQU&-EQUU,'T 1 LA (91 D 82)! 


SALESMODEL4*DMODFL 'EQU1,FQU2,EQU3,FQU4' 


“ትይ 


214+(227,177 T W1 D W2) : (.546535 T 1 LAG S1) P (.1202u8 T DI) 
$2-( 87.0728 E: W1 D #2) P (.829077 T 1 LAG S2) P (.121566 T DI) 


W1+(.304443 T 1) P (.333374 T S1 D 52) P ( .188266 T 1 LAG S1 D S2) 
W2+(.385065 T 1) P (7.291008 T S1 D S2) P (.198681 T 1 LAG #1 D 32) 


A SIMULATION HORIZON 


11958 1 ST 1 1965 1 e 


SSOLVE SALESMODEL4 


S1 


195901 1377.148742 1528.039927 1560. 921662 1594.867777 1617.631134 
1715.438371 1770.080409 1702.184749 


W1 


196001 0.4535188862 0.4547082702 0.4602891791 0.4738244781 
0. 4944205769 0.5129434986 0.4960040896 


)COPY EPLDEMO BASSDATA 
SAVED 9.31.06 01/22/76 


SSOLVEF SALESMODFLY 


91 
195901 1377.148742 1507.507012 1550.708023 1583.611901 1601. 253612 
1684.704856 1744.488829 1741.689299 

W1 


195901 0.451 0.4535188862 0.4547082702 0.4602891791 0.473824478141 
0. 4944205769 0.5129434986 0.4960040896 @ 
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5. DATA FILE ROUTINES 


| This section describes file usage under APLSV. It also 
applies to file usage under APL/CMS with the changes noted in 


5.4 below. 
5.1. SETUP OF A FILE 


The APL ECONOMETRIC PLANNING LANGUAGE provides the 
user with functtons for storing, updating and retrieving 
data Into or from a disk storage device. Data can be either 
time series, cross-sectional series, equations or models. 
The file functions do not distinguish between form of 
representation or rank of the data, but consider each data 
type uniformly as one object. 


A file with an optional NAME ts set up by the function 


DDATA 'NAME' 


The function asks for a maximum number of objects to be 
stored and a "maximum number of elements" per object. The 
second item is used in the allocation of disk space for each 
object. It should be specified in the following manner: 

l. If it is assumed that massive demand for space will arise 
primarily from long time series, the number of values for 
the largest series should be supplied. 

2. If equations or models are relatively larger, 8 times the 
number of characters needed to represent the largest object 
should be supplied. However, it is usually not economical to 
store a multi-equations model as only one object, but it is 
preferable to store the component equatlons separately. 
After retrieval, the model can be reassembled with the 
function DMODEL; see 2.4. 


As with all other functions in this section, error 
messages refer to TSIO messages. For their encoding the user 
should consult the APLSV User's Guide (SH20-1460). 


5.2. DATA STORAGE, UPDATE AND RETRIEVAL 


- Once a-data file is set up with a specific NAME, 
objects can be stored by means of the function 


'NAME' STORE LIST 
LIST contalns the names of the objects, separated by commas. 
= Objects in LIST which were already stored on the file are 
updated, i.e. the previous contents are replaced by the new. 
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A message LISTA NOT SUCCESSFULLY STORED is printed If 
the storing of part of the objects, identified in LISTA, 
falls. In addition various TSIO error conditions may occur. 
If the storage failed due to space allocatton problems a 
new, additional file should be set up. A WN NOT DEFINED 
message indicates that an object NN was not found in the 
workspace, 


Retrievals from a file NAME are obtained by 


'NAME' RETRIEVE LIST 


LIST contains the names of the objects to be retrieved, 


An EMPTY DATASET ERR occurs if retrieval ts attempted 
from a data file which is set up but contains no objects. A 
NN NOT FOUND message ts delivered if the object WN was not 
found tn the file. 


Example 20 


5.3. AUXILIARY FILE FUNCTIONS 


The function 
*WAMENEW ፣ RDATA 'NAMEOLD' 


Men the user to change the name of an existing data 
e. 


The function 


CDATA 'NAME' 


returns the names of all objects stored in NAME. An EMPTY 
DATASET ERR indicates a set up file with no loaded object. 


To erase a data file NAME, one executes 


EDATA 'NAME' 





5.4 FILE ROUTINES UNDER APL/CMS 


File usage under APLSV applies also to APL/CMS, but with 
the following changes: 


The APLSV file routine names have a "C" added under APL/CMS: 
DDATAC, STQREC, RETRIEVEC, RDATAC, CDATAC, AND EDATAC. 
In DDATAC the user is not prompted for size of space needed. 


The error messages are generated by AP110, and are listed in 
the APL/CMS User's Manual (SC20-1846). 
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- DDATA 'SLSMODEL' 
. MAX NO OF OBJECTS: 20 
MAX NO OF ELEMENTS: 25 


n (THE NAME OF THE FILE SHOULD NOT BE LONGER 
a THAN 8 CHARACTERS!) 


LIST1+'E1,E2,E3 .E4,E5,E6! 


LISTI ASSIGN pene eee 


E1 | | i 
S1<(227.177 Z Wi D #2) P P (.546535 T 1 LAG S1) P «120218 T DI 


LISTeLISTi,! ,51,82,1,N2' 
LIST 
E1,E2,E3,Eu ,E5,E6,51,92 ,W1,W2 


"SLSMODEL' STORE LIST 
)ERASE E1 S1 
'SLSMODEL' RETRIEVE 'S1,E1' 


E1 | | | 
S1<(227, 177 T Wi D 929 P X. 546535 T 1 LAG $1) P «120248 T DI 


m——— 
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6. ADDENDUM 


The addendum contains four additfons to the system , 
with the functional objectives: 


1. retranslation of system operators into APL operators, 
2. output of tables with text, 
3, updating of time series, 


ዛ, adjustment of pertodictties for time series. 
6.1. RETRANSLATION: RETRANS 


The function translates operators from the system 
representation to the APL representation. Its argument can 
be the character representation of either one equation or of 
a model with system operator symbols. See example 21. 


The retranslation function calls a function VOCABULARY which 
contains a list of the vocabulary currently used. This 
vocabulary can be altered or extended subject to the 
following rules: 


The APL operators on the right hand side of the 
vocabulary must not be longer than the system operators. 


In case they are shorter, (e.g., PW might be translated 
into x), the APL character must be used as ` 
fill-character to make the lengths of the two operators 
equal. 
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A EXAMPLE 21 


E1+'A+B P C PW (A D B) M (1 LAG A)? 
F2<'C+NQZ A AND B' 


MOD+DMODEL 'E1,E2' 


E1 
A+B P C PW (A D B) M (1 LAG A) 
E2 
C-NOT A AND B 
MOD 
A+B P C PW (A D B) M (1 LAG A) 
C-NOT A AND B 
RETRANS E1 


A+B + C * (A + B) - (1 LAG A) 


RETRANS E2 
C+" A^ B 


RETRANS MOD 
A-B * C * (A * B) - (1 LAG 4) 
C+" A ^B 
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6.2. OUTPUT OF TABLES WITH TEXT: TABULTEXT 


This function is a minor extension of TABULATE. Using 
the same arguments as TABULATE it allows the use of a header 
and the insertion of text in place of variables for which 
row information is to be printed. The information is passed 
by means of global variables defined by the user. 


The desired header is placed into the variable: TABHEAD. 


Row text in place of a variable name is placed into 8 
variable ZTEXT(name). E.g, text for the variable 51 should 
be put into the variable ZEXTS1. 


When the variables  TABHEAD and TEXT... do not exist, the 
output of TABULTEXT is the same as that for TABULATE. 


See example 22 on the next page. 
6.3. UPDATING OF TIME SERIES: UPDATE 


It may be desirable to alter values of time series or 
to insert new values. UPDATE permits: 
. updating of coefficients within the range of the time 
series, | 
. extension of the old time series. 


The use of UPDATE will be made clear by example 25. 
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LABHEAD*'SALES PEK CAPITA ΟΡ FILTER AND NONFILTER CIG' 
LABHEAD-ZABHEAD, "ARELLES BETWEEN 1953 AND 1965! 


TEXISi1-'SALES BRAND 1] 
LEXUS 2+! BRAND 2' 


21 

195401 250.328 398.658 690.081 952.423 1232.111 1414.718 
1526.195 1571.303 1613.578 1657.487 1731.53 
1672.207 1702.119 


S2 
195401 3196.318 2856.406 2668,82ዛ 2453.522 2178.211 


2008.075 2006.827 1991.59 1971.696 1900.684 
1808.007 1590.523 1633.202 


15 0 9 1 LABULYEXY 'S1,S2' 


SALES PER CAPITA OF FILTER AND NONFILTEH CIGARETTES BETWEEN 1953 
AND 1965 


1953 1951 1955 1956 1957 
SALES BRAND 1 250.3 398.7 690.1 952.4 1232.1 
BRAND 2 3196.3 2856.4 2668.8 2453.5 2178.2 
1958 1959 1960 1961 1962 
SALES BRAND 1 1515.7 1526.2 1571.3 1613.6 1657.95 
BRAND 2 2008.1 2006.8 1991.6 197147 1900.7 
1963 1964 1965 
SALES BRAND 1 1731.5 1672.2 1702.1 
BRAND 2 1808.0 1590.5 1633.2 








A EXAMPLE 23 
TT+4 1950 2 DF 18 
TT 

78020012 3556 78. 


DISPLAY 'TT' 


VARIABLE TT 


PERIODICITY = 4 

ORIGIN = 1950 2 

NO OF ENTRIES = 8 

TIME VALUE TIME 

1950 2 1.00000 1951 
1950 3 2.00000 1951 
1950 5 3.00000 1951 
1951 1 4.00000 1952 


| TTT-UPDATE TT 

FIRST PERIOD TO BE UPDATED 
0: | 
u 1951 2 


NEW VALUE(S) 
D: 
1.5 2,5 


RESPECTIVE OLD VALUE(S): 
5 6 


DISPLAY 'TTT' 


VARIABLE TTT 


PERIODICITY = 4 
ORIGIN = 1950 2 
NO OF ENTRIES = 8 


TIME VALUE TIME 
1950 2 1.00000 1951 
1950 3 2.00000 1951 
1950 4 3.00000 1951 
1951 1 4.00000 1952 


P £ Q ኮን 


H|P τω N 
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VALUE 


5.00000 
6.00000 
7.00000 
8.00000 


VALUE 


1.50000 
2.50000 
7.00000 
8.00000 
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TTT+UPDATE TT 
FIRST PERIOD TO BE UPDATED 
U: 
| 4 1949 2 | | 
UPDATING POSITION BEPORE BEGINNING OP SERIES, 
NEW VALUE(S) 
U: | 
1.5 2.5 3.5 4,5 
RESPECTIVE OLD VALUE(S): 
0000 
DISPLAY 'TTT'! 


VARIABLE TTT 


PERIODICITY = LW 

ORIGIN = 1949 2 

NO OF ENTRIES = 12 

TIME VALUE TIME VALUE 

1949 2 1.50000 1950 4 3.00000 
1949 3 2.50000 1951 1 5.00000 
1949 4 3.50000 1951 2 5.00000 
1950 1 5.50000 1951 3 6.00000 
1950 2 1.00000 1951 4 7.00000 
1950 3 2.00000 1952 1 8.00000 


TTT-UPDATE TT 
FIRST PERIOD TO BE UPDATED 
Li: 
4 1949 4 
UPDATING POSITION BEFORE BEGINNING OF SERIES. 
NEW VALUE(S) 
U: 


9 8 7 6 
RESPECTIVE OLD VALUE(S): 
001 2 


DISPLAY 'TTT' 


VARIABLE TTT 


PERIODICITY = 4 

ORIGIN = 1949 4 

NO OF ENTRIES = 10 

LIME VALUE TIME VALUE 

1949 4 9.00000 1951 1 4.00000 
1950 1 8.00000 1951 2 5.00000 
1950 2 7.00000 1951 3 6.00000 
1950 3 6.00000 1951 4 7.00000 
1950 4 2.00000 1952 1 8.00000 
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ITT-UPDATE TY 
"FIRST PEHIOD TO BE UPDATED 
Li: 
4 1952 2 
NEW VALUE(S) 
Lis 
1.6 9.3 
RESPECTIVE OLD VALUE(S): 
0 0 | 


DISPLAY 'TTT' 


VARIABLE TTT 


PERIODICITY = 4 
ORIGIN = 1950 2 
NO UF ENTRIES = 10 
TIME VALUE TIME VALUE 
1950 2 1.00000 1951 3 6.00000 
1950 3 2.00000 1951 4 7.00000 
1950 ዛ 3.00000 1952 1 8, ሀሀሀሀሀ 
1951 1 4,00000 1952 2 ቷ,6ሀሀሀሀ 
TTT+UPDALE T 
FIRST PERIOD TO BE UPDATED 
Lj: 
4 1952 4 
GAP BETWEEN FIRST PERIOD TO BE UPDATED AND LAST PERIOD OF THE ULD 
TIMESERIES. 


INPUT FOR THESE PERIODS IS ASSUMED TO BE ZERO. 
NEW VALUE(S) | 


L: 

11.9 96.3 
RESPECTIVE OLD VALUE(S): 
0 0 

DISPLAY 'TTT' 


VARIABLE TIT 


PERIODICITY = 4 
ORIGIN = 1950 2 
NO OF ENTRIES = 12 
TIME VALUE TIME VALUE 
1950 2 1.00000 1951 4 7.00000 
1950 3 2.00000 1952 1 8.00000 
1950 4 3.00000 1952 -2 . 00000 
1951 -1 4.00000 1952 3 «00000 
1951 2 5.00000 1952 4 11.90000 
1951 3 6.00000 1953 1 96.30000 
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6.4. ADJUSTMENT OF PERIODICITY: PREP 


The right argument of PREP specifies the time series 
to be altered, the left argument gives the new periodicity. 
The adjustment may either correspond to a compression or to 
an expansion (interpolation) of the time series. For each 

time series in the right argument, then, one has to specify 


the desired method of compression or interpolation. 


We associate a variable PERCHAN(name) with the time 
series (name). If this variable exists, it [5 used during 
execution of PREP, if not it is generated through a question 
at the terminal. The first component of PERCEAN(name) 
specifies a method for compression, the second a method for 
interpolation. It will always be clear from the context 
which one is to be used. 


Consult the description of CHANGE for possible error 
messages. Various possibilities are demonstrated in example 
2h. | 
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EXAMPLE 24 


Ti«4 1950 3 DF 1 2 3 ዛ 5 6 7 8 9 

T2+12 1960 7 DE 130 
12 | 

2352712 12345 6 7 8 9 10 11 12 19 14 15 16 17 18 
19 20 21 22 23 24 25 26 27 28 29 30 


2 PREP 'T1,T2' 

STATE COMPRESSION AND INTERPOLATION METHOD FOR T1 
3 2 

PERCHANT1: 3 2 

STATE COMPRESSION AND INTERPOLATION METHOD FOR T2 
21 ~~ 

PERCHANT2: 2 1 


T1 
390202 24 6 8 9 
T2 
392202 21 57 93 129 165 
EC T1[1] 
2 1950 2 
EC T2[1] 
2 1960 2 


12 PREP 'T1,T2' 
PERCHANT1: 3 2 
PERCHANT2: 2 1 
T1 
2340712 0.3333333333 0.6666666667 1 1.333333333 1.666666667 
2 2.333333333 2.666666667 3 3.333333333 3.666666667 
4 4.333333333 4.666666667 5 5.333333333 5.666666667 
6 6.333333333 6.666666667 7 7.333333333 7.666666667 
8 8.166666667 8.333333333 8.5 8.666666667 8.833333333 
9 
T2 | 
2352712 21 21 21 21 21 21 57 57 57 57 57 57 93 99 
93 93 93 93 129 129 129 129 129 129 165 165 
165 165 165 165 
EC T1[1] 
12 1950 7 
EC T2[1] 
12 1960 7 


52 PREP 'T1,T2' 
PERCHANT1: 3 2 
CHANGE PER ERR 
WITH SERIES: T1 
PERCHANT2: 2 1 
CHANGE PER ERR 
WITH SERIES: T2 
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